Andreev Level Qubit 



A. Zazunov a , V. S. Shumeiko a , E. N. Bratus' 6 , J. Lantz a , and G. Wendin Q 
a Department of Microelectronics and Nanoscience, Chalmers University of Technology and Goteborg University, 

S-41296 Goteborg, Sweden 
b B. Verkin Institute for Low Temperature Physics and Engineering 
310164 Kharkov, Ukraine 

We investigate the dynamics of a two-level Andreev bound state system in a transmissive quantum 
point contact embedded in an rf-SQUID. Coherent coupling of the Andreev levels to the circulating 
supercurrent allows manipulation and read out of the level states. The two-level Hamiltonian for 
the Andreev levels is derived, and the effect of interaction with the quantum fluctuations of the 
induced flux is studied. We also consider an inductive coupling of qubits, and discuss the relevant 
SQUID parameters for qubit operation and read out. 
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Recent observations of quantum coherence in super- 
conducting circuits 0-^) have made superconducting 
qubits a realistic possibility. Superconducting qubits em- 
ploy the phenomenon of macroscopic quantum coherence 
(MQC) 0, and operate with coherent superpositions of 
quantum states of a macroscopic object - a superconduct- 
ing condensate. An elementary MQC circuit consists of 
a hysteretic SQUID with a small capacitance Josephson 
tunnel junction (persistent current or flux qubit). The 
qubit is biased at half-integer external flux, and oper- 
ates with the two degenerate current states corresponding 
to the clockwise and counter-clockwise circulating persis- 
tent currents. The coupling of the current states is due to 
quantum fluctuations of the electric charge on the junc- 
tion capacitor. Quantum measurements of the fluctuat- 
ing persistent current, or induced magnetic flux, provide 
means to read out the qubit state. Further modifica- 
tions of the flux qubits involve implementation of mul- 
tiple junction circuits In multiple junction MQC 
circuits, it is possible to employ one of the available dy- 
namic variables for qubit operation, and another one for 
qubit readout 

In this paper we consider a new type of superconduct- 
ing qubit where the switching between the two persistent 
current states in a SQUID is achieved by employing a 
true microscopic system formed by the two-level Andreev 
bound states in a superconducting atomic-size quantum 
point contact embedded in the SQUID |)|. In this An- 
dreev level qubit, the quantum information is stored in 
the microscopic quantum system, the Andreev bound 
states, similar to non-superconducting solid state qubits 
like localized spins on impurities ]To[ | or quantum dots 
jn]]. Read-out of the Andreev level qubit is achieved 
by monitoring the macroscopic persistent current or the 
induced flux in the SQUID, similar to the MQC qubits. 

The Josephson effect in a single atomic-size quantum 
point contact (QPC) embedded in a low-inductance non- 
hysteretic SQUID (Fig. 1) has first been investigated 
by Koops et al. Gaj; in this experiment, the averaged 



current-phase relation in the ground state was measured 
by performing a classical measurement of the induced 
flux. The results of this experiment, and also of other 
experiments on atomic-size QPCs where the critical cur- 
rent [jl3| and current- voltage characteristics |l4|jr^ ] have 
been investigated, are consistent with a theoretical pic- 
ture in which the Andreev bound levels play the central 
role in the Josephson current transport DM . 




FIG. 1. Left: Sketch of the Andreev level qubit - a 
non-hysteretic rf-SQUID with a quantum point contact 
(QPC), and the equivalent circuit containing the Josephson 
junction and LC-oscillator. Right: The energy spectrum of 
the QPC with finite reflectivity (R — 0.04) (solid line), ap- 
pears as a hybridization of the current states (R = 0) (dashed 
line) . 

The Andreev bound levels are formed in a QPC due 
to Andreev reflections by the discontinuity of the super- 
conducting phase at the contact in the presence of the 
applied current. The Andreev level wave functions are 
localized in the vicinity of the contact over a distance of 
the order of the superconducting coherence length, and 
the number of Andreev bound levels is limited to one pair 
of levels per conducting electronic mode. Thus, a super- 
conducting QPC may be viewed as a kind of quantum 
dot which contains a finite number of localized quantum 
states. In highly transmissive QPCs, the Andreev levels 
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lie deep within the superconducting gap and are well de- 
coupled from the continuum quasiparticle states in the 
electrodes. 

The most important property for qubit applications is 
the coherent coupling of the Andreev levels to the super- 
current flowing through the contact. This makes the An- 
dreev levels accessible for manipulations and for measure- 
ments. Extensive studies have shown that the time evolu- 
tion of the Andreev levels can be controlled by applying 
resonant rf flux pulses JTij ] or by ramping the external 
flux P). By these means, one can drive Andreev levels 
out of the ground state and prepare any excited state 
|l8l| . On the other hand, the time evolution of the An- 
dreev levels changes the current through the QPC. Thus 
the QPC operates as a quantum switch which controls 
the direction of the circulating current in the SQUID. 
Performing quantum measurement of the circulating cur- 
rent, or of the corresponding induced flux, in the SQUID 
one is able to measure the state of the Andreev level sys- 
tem. Fidelity of such a measurement requires the qubit 
evolution to be slow on a time scale of intrinsic electro- 
magnetic fluctuations in the SQUID. Then the current 
(and induced flux) averaged over electromagnetic fluctu- 
ations will adiabatically follow the qubit evolution. 

It is important to note that the interaction of the An- 
dreev levels with electromagnetic modes in the SQUID 
is an essential element of the qubit dynamics: To main- 
tain the current switching, the plasma frequency of the 
SQUID must be sufficiently large. Furtheremore, since 
the current in the single-mode QPC undergoes strong 
quantum fluctuations, the induced flux and hence the 
superconducting phase difference, are also fluctuating 
quantities, and the theory of the Andreev level qubit 
should include a full quantum mechanical treatment of 
the coupled Andreev levels and electromagnetic fluctua- 
tions. 

To derive an effective quantum Hamiltonian describing 
coupled Andreev levels and electromagnetic fluctuations, 
we employ a path integral approach commonly used in 
MQC theory The central problem here is to extend 
the theory, originally developed for tunnel junctions, to 
the for us interesting case of high transmission QPCs; 
this problem is solved by using the exact boundary con- 
dition in the action instead of the tunnel Hamiltonian. 
Following Ref. |l9] , we present the evolution operator for 
the system on the form, 



U = Vcj)e 
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where the integration is performed over the super- 
conducting phase difference 4>(t) and Nambu-Grassman 
fields 4 ! L,R(r,t) representing electronic degrees of free- 
dom in the left and right electrodes respectively. The 
action S osc [0] describes an LC-oscillator formed by the 
junction capacitance C and the superconducting loop in- 
ductance L (see Fig. 1), while the action Sj describes 



the Josephson junction. 

To derive the action for the Josephson junction, we 
consider the two superconducting reservoirs coupled via a 
single-channel localized scattcrcr. The scatterer is repre- 
sented by a normal electron scattering matrix, which im- 
poses the boundary condition for the quasiclassical wave 
functions of the quasiparticles in the reservoirs. The mi- 
croscopic Hamiltonian for the reservoirs, 

H= J2 [dr^M (T ,h=(^--f?ja z +Aa x , (2) 

is considered within the mean field approximation, the 
superconducting phase being gauged out and included 
in the boundary condition. The Hamiltonian gener- 
ates the Lagrangian, C — J2a=L u. I ^ r ^a^fo^t — h)ip<r- 
In the boundary condition, which connects the fields 
ipL,n(0,t) at the contact, the energy dependence of the 
scattering matrix on the scale of A can be neglected, 
so that the same matrix describes both the electrons 
and holes. Furthermore, without loss of generality, it 
is possible to eliminate constant scattering phases from 
the boundary condition and include them in the posi- 
tive and negative momentum components of the fields 
?/v(r, i). Within the quasiclassical approximation, such 
a transformation will not affect the Hamiltonian in Eq. 
(||) . The boundary condition can then be written on the 
form, X (t) = de-^Wy^L&t) - (I + r)if) R (0,t) = 0, 
where d = \f~D and r — \J~R are real transmission and 
reflection amplitudes respectively. This boundary condi- 
tion is incorporated into the junction action Sj by means 
of the Lagrange field ?/(i), 



JSj/h 



V 2 n e 



(i/h) J dt(c-nx-xv) 



(3) 



We will sec later that the fcrmionic field rj (t) describes 
the Andreev levels. 

Integration over the rapidly varying fcrmionic fields in 
the electrodes, tp a (r,t), yields the effective action 



iS ef[ [<t,]/n _ f V 2 (i/h) j dt t dt 2 r)(ti)G(ti,t 2 h(t2) 
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-<?(*! -* 2 ), g(u) = J2(^-h)- 1 - (5) 



Equation (^) provides the required generalization of the 
the effective action of the tunnel theory |L9| to junctions 
with arbitrary transparency. In the low frequency limit, 
uj <C A/?L, the Green's function g(to) reduces to a simple 
form, g(tu) = (—Tti>F/A)(h~ui + Ao- x ) (yp is the density of 
states at the Fermi level), which leads to the local-in-time 
effective Lagrangian, G{t\,t2) = G(ti)S(ti — t 2 ), 
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G(t) = ihdt 



hd, 



U ze 



-i0/2 



(6) 



with z = A[cos(0/2) + ir sin(0/2)]. 

Equation (^|) describes a two-level fermionic sys- 
tem. Under stationary conditions, d t cf> = 0, the 
spectrum of the system, huj = ±E a , E a — 



Ay cos 2 (0/2) + i? sin 2 (0/2), coincides with the Andreev 
level spectrum jL6|, shown in Fig. 1. It follows from this 
equation that the assumed low-frequency approximation 
is appropriate for transparent contacts, R -C 1, at (f> rs tt, 
where the Andreev level energy is small, E a <C A. To de- 
rive the Hamiltonian for the Andreev two-level system, 
it is convenient first to eliminate the time derivative of 
the phase in Eq. (|J) by means of a unitary rotation, and 
then to use the relation, G(t) = ihdt — H a , to obtain 



H a = 







_ z * £ ir4>/2 



-ze~ ir 't'/ 2 




(7) 



This equation is the first main result of the present paper 

The slow dynamics in transparent contacts, R <C 1, 
is described by two variables, 77(f) and <p(t). This is 
rather different from tunnel contacts with D <C 1, where 
the effective action only depends on the phase difference 
and basically reduces to the potential Josephson energy 
JL9[ . This difference is easily understood if one takes into 
account that the Andreev levels in tunnel junctions are 
close to the gap edge, E a w A, making r/(t) a rapid vari- 
able. Integration over 7/ in Eq. (^]) assuming small D in 
Eq. (||) , recovers the effective action of Ref. [[L9| . 

Let us now consider the current through the junction. 
The statistically averaged current (I) can be expressed 
through the Josephson part, Uj, of the evolution opera- 
tor in Eq. (g), (I) = 2ei(5/5cf>)hx(Uj). In terms of the 
effective action, the equation for the averaged current re- 
duces to the form (I) = Tr{p rt I), where p v is the density 
matrix, and / is the current operator of the two-level 
Andreev system, 



2e dH a el{(f>) 



e 
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T{4>) = AD sm(<f>/2). The current operator I does not 
commute with the Hamiltonian H a , which is a conse- 
quence of the normal electron reflection at the QPC. 
Therefore the Andreev levels consist of superpositions of 
the current eigenstates, unless R = 0, and hence the 
current expectation value in the Andreev state, I a = 
(2e/h)(dE a /dc/)) = (e/2h){DA 2 / E a ) sin0, does not co- 
incide with the current eigenvalues, ±eT/h, which are 
evaluated during the quantum measurements. Further- 
more, the Andreev level current undergoes quantum 
fluctuations with the spectral function (cf. Ref. |H|), 



S(u) = I 2 a Rt&n 2 {4>/2)5{uj - 2E a ). In the SQUID geome- 
try, these fluctuations generate strong quantum fluctua- 
tions of the phase. 

We now take the quantum dynamics of the super- 
conducting phase into the consideration. The quan- 
tum Hamiltonian of the iC-oscillator associated with 
the action S asc [<f>\ in Eq. ([!]) has the form, H osc — 
-(hd 4> ) 2 /2M + Mtu 2 j> 2 /2, where M = h 2 /8E c , w = 



8E L E c /h 2 , E c = e 2 /2C, E L = {hc/2e) 2 (l/ L), and 

(j> is related to the induced flux, $ = (ch/2e)4>. 

In the practically important case of small loop induc- 
tance, Ej -c El, the SQUID is in the non-hysteretic 
regime, and the induced flux is small, <f> <C 1. Introducing 
the phase difference <p e related to a stationary external 
flux we expand the Hamiltonian in Eq. (j?]) over small </>; 
then the Hamiltonian of the whole system in the current 
eigenbasis takes the form, 



H = - 



A I < ( >s ^ a z + r sin — a. 



H 



(9) 



This Hamiltonian describes a spin degree of freedom lin- 
early coupled to an oscillator, the steady state of the 
oscillator being shifted from the origin by ±T/2Mui 2 de- 
pending on the spin direction (direction of the current 
in the junction). We are interested in the case when the 
induced flux adiabatically follows the evolution of the 
Andreev levels. This regime corresponds to a large oscil- 
lator frequency compared to the Andreev level spacing, 
hui ^> 2E a . In this case, the oscillator can be assumed to 
be in the ground state, <fio(4>±), <fi± — (f>±T/2Muj 2 , since 
the probability of transitions among the oscillator levels 
is small. Averaging out the ground state phase fluctua- 
tions P^l , we finally arrive at the effective Hamiltonian 
describing the Andreev level qubit, 



H q = — A ( cos a z + qc>r sin a x 



(10) 



The factor, go = exp (—X 2 / iMhui 3 ) , is the overlap in- 
tegral between the oscillator ground state wave func- 
tions for different current directions, ipo(cj>±). The av- 
eraged value of the induced phase is then given by 
(4>) = (l/2Mu! 2 )Tr(p v a z ), which yields the relation be- 
tween the induced flux operator and the current opera- 
tor, $ = (L/c)I. Therefore quantum measurement of the 
flux on a time scale l/ui < r < h/2E a will allow correct 
evaluation of the Andreev level state. 

The qubit Hamiltonian ( |Io| ) is another main result of 
this paper. Equation (10) is equivalent to the Hamilto- 
nian of non-interacting Andreev levels, Eq. (j^), but with 
reduced reflectivity, R = QqR. One may interpret this re- 
duction as the effect of the inertia of the loop oscillator, 
which makes it more difficult for the Andreev levels to 
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switch direction of the current. The effect becomes in- 
creasingly strong in the limit of a classical oscillator with 
large " mass" . This renomalization effect leads to reduc- 
tion of the Andreev level energy E a — > E a , and hence 
to the reduction of the frequency of the qubit rotation, 
remaning the results of Ref. Jl^Jl^l essentially the same. 
This might be important for practical applications, be- 
cause it would allow the frequency of the qubit rotation 
to be tuned by choosing circuit parameters rather than 
by tuning the contact reflectivity. 

The following chain of inequalities summarizes the re- 
quirements for the Andreev level qubits, 

2Av / l'< hu> < A E L . (11) 

These requirements are consistent with the typical cir- 
cuit parameters of the experimental MQC qubits HQ. 
The critical current in the QPC is close to the maxi- 
mum supercurrent of a single mode, I c ~ eA/h (~ 400nA 
for Nb), and the Josephson energy, Ej A. Assuming 
L ^O.lnH, and C ~ O.lpF, we estimate u> ~ 10 11 sec _1 , 
and El/% ~ 10 13 sec _1 ~ 10 A^/?),. For these values, the 
<7o-factor is of order unity, and the reflectivity of the con- 
tact must be small, R < 0.01, which is, in principle, ac- 
cessible in the experiments with atomic-size QPC |l3| , |l5| . 
However, even a slight increase of the inductance will sig- 
nificantly decrease the go-factor, and the constraint on 
the bare contact reflectivity will become less restrictive. 
In summary, we estimate the upper bound for the qubit 
operation frequency to be about 10 10 sec _1 . 

The qubit operation frequency must significantly ex- 
ceed the relaxation and dephasing rates of the qubit. 
The relaxation and dephasing mechanisms (external flux 
fluctuations, radiation, etc.), which have been extensively 
discussed for the MQC flux qubits Q , are also relevant for 
the Andreev level qubit. The interaction of the Andreev 
levels with microscopic degrees of freedom in the contact, 
primarily with the phonons, does not impose any further 
limitations. Investigation of this problem has shown ]23j ] 
that the relaxation rate is very sensitive to the Andreev 
level spacing: for small spacing, t _1 ~ R 2 r~ h (A) at 

small temperature, T <C E a . For R < 0.01 this relax- 
ation rate is smaller than the qubit operation frequency 
by at least a factor of 10 4 . 

We conclude with a discussion of the qubit-qubit in- 
teraction. Let us consider an inductive coupling as the 
most relevant interaction for the flux qubits. This cou- 
pling will introduce hybridization of the loop oscillators, 
which is described by inserting the inductance matrix in 
the H osc term in Eq. (^|). The averaging over flux fluctu- 
ations is conveniently performed using the oscillators nor- 
mal modes. The averaging procedure leads to the qubit 
Hamiltonians in Eq. (|l0|) with slightly different dress- 
ing factors, and to a Hamiltonian of direct qubit-qubit 
interaction. For the two qubits the effective interaction 
has the form, Hi nt — —{ejcK) 1 ( y MX{l2) o z \o Z 2, where 



M. is the mutual inductance. The two-qubit configura- 
tion may consist, in particular, of a single QPC with two 
conducting modes; in this case M. = L. 

In conclusion, we have developed a theory for the An- 
dreev level qubit, a device consisting of a SQUID with a 
quantum point contact, combining the features of mi- 
croscopic and macroscopic quantum systems. We de- 
rived the two-level Hamiltonian for the Andreev levels 
and showed that it is strongly dressed by the quantum 
fluctuations of the induced flux. We also derived the ef- 
fective interaction Hamiltonian for inductively coupled 
qubits, and discussed the relevant circuit parameters for 
the qubit operation and read out. 
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